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I present predictions for the inclusive production of D mesons at the CERN LHC in the 
general-mass variable-flavor-number scheme at next-to-leading order. Numerical results 
are compared to data where available, and uncertainties due to scale variations, parton 
distribution functions and charm mass are discussed. I point out that measurements at 
large rapidity have the potential to pin down models of intrinsic charm. 

D-mcson production at the LHC was studied by the ALICE Q], ATLAS 0, and LHCb 
Collaborations [3j ■ Here I present predictions for the inclusive production of D mesons at the 
LHC within the general-mass variable-flavor- number scheme (GM-VFNS) [4]. More results and 
additional details of the calculation can be found in [5]. In a recent paper [6J, we have also 
considered the inclusive production of B mesons, for which experimental results from the CMS 
Collaboration are available [7]. For an alternative approach, see Ref. [5]. 

Figure [l] shows a comparison of the GM-VFNS 
predictions for the transverse momentum distribu- 
tion with data from ALICE. Here the renormaliza- 
tion (fin) and factorization scales for initial state 
(fii) and final state (pf) singularities are fixed by 
/ii = £,i\/p T + ra 2 c , where m c is the charm quark 
mass, and the scale parameters & (i = R, F, I) 
are varied about the default values of 1 by factors 
of 2 up and down to obtain an estimate of a the- 
ory uncertainty band (dotted lines in the figure). 
The data are reasonably well described by theory 
at the larger values of px, where data are avail- 
able, but due to the choice = 1, theory starts 
to overshoot at pt < 5 GeV. There, the fixed fla- 
vor number scheme (FFNS) [5] works better (see 
the dashed lines in Fig. [I). The GM-VFNS is pre- 
ferred at large px since it includes resummed con- 
tributions from large logarithms by virtue of the 
DGLAP evolution equations for the parton distri- 
bution (PDFs) and fragmentation functions (FFs). 
The GM-VFNS also predicts smaller scale uncer- 
tainties than the FFNS. We have used CTEQ6.6 
PDFs [15] and, in the case of the GM-VFNS, FFs 
of Ref. [10]. The FFNS calculation is performed 
without including a FF; the transition from the 
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Figure 1: da/dpx for p + p — >• D° + X in- 
tegrated over rapidity in the range — 0.5 < 
y < 0.5 for Vs = 7 TcV at NLO in the GM- 
VFNS (solid line) and the FFNS (dashed 
line) . Dotted lines describe the correspond- 
ing error bands from scale variations. The 
ALICE data were taken from p] . 
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Figure 2: Ratios of da/dpx for D-meson production at ALICE at yfs = 7 TeV using = £f = 
0.8 and £r = 1. All cross sections and the data from Ref. [T] are normalized to the GM-VFNS 
prediction with & = 1. The PDFs are taken from MSTW08-NLO [T3] and the charm quark 
mass is m c = 1.5 GeV. 



charm quark to the charmed meson is taken into account by multiplying the parton level result 
with the branching ratio BR(c — > D°) = 0.628. 

The uncertainties due to variations of the factorization scales are dominant. It is interest- 
ing to see that the scale parameters can be chosen to bring the GM-VFNS predictions into 
agreement with the data also at low values of px- This is shown in Fig. [2] for MSTW08-NLO 
PDFs [13] and using m c = 1.5 GeV for the charm quark mass. The differential cross sections 
da I dpx are shown here for £/ = £p = 0.8, £r = 1 in pt bins and compared with data points 
from the ALICE collaboration pQ . All results are normalized to the GM-VFNS prediction with 
Ci — 1. One can see that a proper choice of the factorization scales can help to ensure that the 
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Figure 3: da/dpr for D- meson production at ALICE at y/s = 7 TeV for different PDFs. All 
cross sections are calculated with £j = £p = 0.7, £r = 1 and normalized to the GM-VFNS 
prediction with CTEQ6.6 PDFs. The histograms from top down correspond to CT10 [llj . 
HERAPDF 1.5 (NLO) Q2], MSTW08-NLO [13] and NNPDF 2.1 [H]. 
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resummed contributions due to incoming heavy quarks, and those due to gluon fragmentation, 
fade out in a controlled manner as pr/m 0, i.e. in the kinematic region where the FFNS 
should be appropriate. 

In Figure [3] an attempt is made to show the uncer- 
tainties coming from using different PDF input. The 
results for most of the bins lie within the error bars 
of the experimental data and do not prefer one PDF 
set over another. Actually, due to different values of 
m c used in the PDF fits and the corresponding dif- 
ferent lengths in the evolution path from the charm 
production threshold up to there is some residual 
m c dependence of the predicted cross sections at low 
values oipx- The value m c = 1.5 GeV used in our cal- 
culation agrees with the one in the fragmentation func- 
tions of Ref. [TO], but not with the one in the parton 
distribution functions used here. While the CTEQ6.6 
and CT10 sets use m c = 1.3 GeV, in the MSTW08- 
NLO, NNPDF 2.1, and HERAPDF 1.5 (NLO) sets 
m c = 1.4 GeV was chosen. A consistent calculation 
would require the same value of m c in all components 
of the cross section formula. However, separate fits of 
the fragmentation functions for different values of m c 
are not available. The dependence on the heavy quark 
mass is, however, not very strong and non- negligible 
only in the low pr range, see Fig. [4] 

Non-perturbative contributions to the charm quark 
content of the proton may lead to enhanced charm par- 
ton distributions c(x,/if) at x > 0.1. This can become 
visible in the cross section for D meson production 
at large rapidities. Parametrizations of this so-called intrinsic charm are available from the 
CTEQ collaboration, based on various models and compatible with the global data samples. In 
Ref. |17] , we have studied the impact of these models on possible measurements at the Tevatron 
and at BNL RHIC. Here, I present results of a calculation using the parametrization CTEQ6.6 
|15) to obtain an estimate of the expected relative enhancements of the pt distributions in bins 
of rapidity. Figure [5] shows typical results for D° production; for other D mesons, the results 
are very similar. Two models have been selected among the possible options in CTEQ6.6 (see 
Ref. [TS] for details): Fig. [5^, shows the calculation using the BHPS model with a 3.5% (c + c) 
content in the proton (at the scale /if = 1-3 GeV), Fig. [5)3 refers to the model of a high strength 
sea-like charm component. In both cases, one observes large enhancements, increasing with ra- 
pidity, and in the first model also with px- Thus one can expect that forthcoming data from 
the LHCb experiment should be able to exclude or narrow down models for intrinsic charm. 
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Figure 4: da/dpx for _D°-meson pro- 
duction at ALICE [I] with MSTW08- 
NLO PDFs [13] normalized to the 
GM-VFNS prediction with £ IiFiR = 1 
and m c = 1.5 GeV. The lower two his- 
tograms are obtained using = 0.7 
and = 1. For the two dashed 
histograms, m c — 1.4 GeV was used 
in both the hard scattering matrix 
elements and in the MSTW08-NLO 
PDFs. 
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Figure 5: Ratio of the p T distributions da/dp T for p + p -> D° + X at NLO in the GM-VFNS 
at s/s = 7 TeV, using two different models of intrinsic charm: (a) BHPS model with 3.5% 
(c+c) content (at fip = 1-3 GeV), (b) model with a high strength sea-like charm component. 
The FFs are taken from Ref. [TU]. The various lines represent the default predictions for 
£h = £j = £f — li integrated over the rapidity regions indicated in the figures (larger rapidities 
correspond to larger cross section ratios everywhere in (a) and at small pt in (b)). 
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